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Abstract 

Using physical insights and advanced first-principles calculations, we suggest that corundum (a- 
AI2O3) is an ideal gate dielectric material for graphene transistors. Clean interface exists between 
graphene and Al-terminated (or hydroxylated) AI2O3 and the valence band offsets for these sys- 
tems are large enough to create injection barrier. Remarkably, a band gap of ~ 180 meV can be 
induced in graphene layer adsorbed on Al-terminated surface with an electron effective mass of ~ 
8 x 10~ 3 m e , which could realize large ON/OFF ratio and high carrier mobility in graphene tran- 
sistors without additional band gap engineering and significant reduction of transport properties. 
Moreover, the band gaps of graphene/A^Oa system could be tuned by an external electric field 
for practical applications. 

PACS numbers: 73.22.-f, 68.65.Pq, 68.43.Bc, 68.47.Gh 







Graphene has been studied intensively due to its unique electronic and mechanical proper- 
ties such as extremely high carrier mobility (over 200000 cm 2 /V-s for suspended samples) 
However, to utilize graphene, which has a zero band gap when it is pure, for electronic devices 
such as the field-effect transistors (FETs), it is essential to open up a band gap in graphene to 
realize the ON/OFF switch function. One way to open the band gap in graphene is utilizing 
the quantum confinement effect, e.g., etching graphene into one-dimensional nanoribbonsj^- 
[5]. In practice, very narrow graphene nanoribbons (GNRs) (~ 10 nm) are necessary to 
achieve a band gap of ~ 0.2 eV (ON/OFF ratio > 10 2 ) [2-5], but a large scale production of 
such narrow GNRs is still quite challenging. Moreover, the carrier mobility of GNRs (about 
hundreds of cm 2 /V-sj4, 5|) is several magnitudes lower than that of graphene sheet due to 
the (intrinsic) band folding and phonon scattering^] as well as (extrinsic) difficulty in con- 
trolling the edge quality in experiments . Another way to open a band gap in graphene 
is breaking the inversion symmetry of the A, B sublattices, e.g., by placing graphene onto 
some special substrate. In this case, the band structure near the Dirac point or the carrier 
mobility of the graphene is better preserved if the interaction between graphene and sub- 
strates is weak. Obviously, this approach has technological advantages over the etching of 
graphene. However, a simple guideline on how to search an ideal substrate that could induce 
a sufficiently large band gap in graphene is still unclear, especially for substrate which can 
be integrated directly into the current FET technology. 

Most of the graphene FETs to date employ silicon oxide (SiOg) as the bottom-gate dielec- 
tric and an alternative high-dielectric-constant (high-AO material as the top-gate dielectric 
(e.g., the amorphous structures of HK^pj] and Al 2 03[8[). The integration of a high-A; top- 
gate can push the FET performance to a much higher limit because it can better screen 
charged impurities and enhance carrier mobility in graphene 7H9[. Analogy to silicon FET, 
an optimal high- k gate dielectric material should have high dielectric constant, large in- 
jection barrier (i.e., large band offset (> 1 eV) with respect to graphene), high chemical 
stability, and no or minimal interface states at the high-&-oxide/graphene interface [10]. In 
addition, an ideal high- k dielectric is desired to have the ability to open an adequate gap in 
graphene for FET operation. 

Although several calculations have shown that a small gap could be induced in graphene 



by Si02 
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12|, Hf0 2 [l3[, or BnQ substrates, a general rule to find a substrate which 



can induce a larger band gap in graphene is still lacked. Here, we propose that one should 
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search for the substrate which has atoms with large chemical potential difference at the 
surface and is commensurate with the graphene lattice. These two conditions can cause large 
potential difference Aab at the A, B sublattices of graphene, thus produce large band gap 
in graphene. Under this physical insight, we suggest that the (reconstructed) Al-terminated 
a-Al 2 03 (0001) surface is an ideal high-k-oxide substrate for graphene FETs. This is because 
after the reconstruction both Al and O are present at the surface and the chemical potential 
difference between Al and O is much larger than that of other popular substrates such as 
SiC>2, BN, and HfC>2 studied before, and the lattice mismatch between CK-AI2O3 (0001) and 
graphene is also relatively small. Carried out by advanced first-principles calculations, we 
demonstrate that a large band gap of ~ 180 meV at the Dirac point appears in graphene 
layer adsorbed on Al-terminated (0001) surface with a quite small effective mass of ~ 8 
x 10~ 3 m e for Dirac fermions. As we expected, the band gap of graphene on AI2O3 is 
significantly larger than on other widely used substrates 15j . The size of the band gaps of 
graphene/ AI2O3 systems could be further tuned by an external electric field for practical 
applications. Moreover, interface states can be eliminated if graphene is grown on Al- 
terminated (or hydroxylated) AI2O3 (0001) surface and the valence band offsets for these 
hybrid systems are large enough for injection barrier. 

All the density-functional-theory (DFT) calculations are performed by using the VASP 
code[16j. Projector augmented wave (PAW) potentials are used to describe the core elec- 
trons, and generalized gradient approximation (GGA) with the PBE functional is selected 
in our calculations. We find that the van de Waals (vdW) interaction plays an indispens- 
able role in accurately determining the adsorption configuration and binding strength in 
this system. The effect of vdW interactions is taken into account by using the empirical 
correction scheme of Grimme (DFT + D/PBE) 17|, which has been proved to be successful 
in describing the geometries of graphene related structuresQ- Our DFT+D/PBE calcu- 
lations show that the interlayer distance of graphite is 3.22 A and the interlayer binding 
energy is -54 meV/C-atom, in excellent agreement with quantum Monte Carlo calculation 
(-56 meV/C-atom) 19[ and experimental value (-52±5 meV/C-atom) 20|]. It is well-known 
that GGA type calculations usually underestimate the band gaps of semiconductors and the 
absolute band edge energy from the GGA calculation is not always reliable. Since hybrid 
functional calculations could give improved results for both conventional semiconductors j2l| 



and graphene nanostructures 



22|, we adopt HSE hybrid functional to calculate the electronic 
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FIG. 1: Side and top views of optimized structures of graphene on (a) Al-terminated, (b) ex- 
terminated, and (c) fully hydroxylated AI2O3 (0001) surface. The green, red, and pink balls denote 
the C, O, and Al atoms, respectively. The yellow lines represent the unit cell in calculations. 

structures of these A1 2 3 /graphene systems. The A1 2 3 (0001) surface is modeled by a slab 
containing six oxygen O3 layers and twelve or eleven aluminum layers (depending on the 
specific surface studied) . The second surface of the slab is passivated by pseudo H atoms and 
a 15 A- vacuum region is included. Single- layer or bilayer graphene with 2x2 periodicity 
is placed on 1 x 1 cell of AI2O3 with a lattice mismatch ~2% 23J]. A T-centered 24 x 24 
x 1 k-point sampling is used for the Brillouin-zone integration. The energy cutoff is set to 
400 eV and structural optimization is carried out on all systems until the residual forces are 
converged to 0.01 eV/A. The dipolar correction has been included 24j. 

a-Al 2 03 has rhombohedral symmetry R3c (D® d , No. 167). The oxygen atoms form a 
nearly hep structure, and the metal atoms fill two-thirds of the octahedral sites between the 
oxygen layers 25j. The (0001) surface of a-Ai20 3 is of major technology importance because 
it is often used as a substrate in growth of semiconducting as well as superconducting 
materials. Low-energy electron diffraction (LEED) revealed that the (0001) crystal surface 



exhibits a 1 x 1 structure below ~1250° in air or in vacuum Although theoretical 

calculations demonstrated that Al-terminated (0001) surface is the most stable one 28|, |29| . 
both Al- and O- terminated (0001) surfaces are observed in experiments and considered in 
our present work 25M27] . 

For the Al-terminated (0001) surface, our calculations show that the topmost Al atoms 
move down ~ 0.65 A into the next oxygen layer after the relaxation compared to the cleaved 



surface, as shown in Fig. la. This surface reconstruction stabilizes the surface by the 
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large suppression of surface polarization [28l. |29|. In order to investigate the stable interface 
structures, we have calculated all the high symmetrical arrangements between graphene and 
Al-terminated surface, where a surface Al or O atom is directly below either: a graphene C 
atom (Tai or To), the hollow site of graphene C atoms (Hai or Ho), or the center of C-C 
bridge site (Bai or Bo). In addition, the case of a C atom directly above the center of surface 
Al-0 bridge site is also considered. Our DFT+D/PBE total energy calculations indicate 
that the most stable configuration is the Tai (Fig. la) with an adsorption energy of -76.4 
meV/C-atom and the distance between graphene layer and substrate surface is ~ 2.74 A. 
The adsorption energy of the T^ configuration is more negative than other high symmetrical 
configurations by 2 ~ 9 meV/C-atom. It should be noticed that the vdW correction is very 
important in this case: without the vdW correction, the calculated interface distance of 
the Tai configuration would be drastically increased to 3.31 A and the adsorption energy 
would be only -2.6 meV/C-atom; the energy difference between Tai and other configurations 
is within 2 meV/C-atom in the absence of vdW correction. LDA calculation could give a 
similar interlayer distance between graphene and substrate surface, but the binding energy 
is still largely underestimated by 34 meV/C-atom. 

On the O-terminated surface, as each O atom has one dangling bond, the O-terminated 
surface is chemically reactive and could strongly interact with the graphene layer [3^]. For 
the lowest energy configuration, all three surface O atoms in the unit cell form chemical 
bonds with the graphene layer to suppress the surface dangling bonds, as shown in Fig. lb. 
The C-0 bond length ranges from 1.45 to 1.50 A. These C-0 binding severely distort the 
planar graphene structure (C-C bond length varies from 1.36 to 1.49 A) and the binding 
energy between graphene and O-terminated surface is -349 meV/C-atom, which is much 
larger than that of Al-terminated surface. In addition, fully hydroxylated AI2O3 surface is 
also considered, as H atoms can be incorporated into the bulk structure during growth 25]. 
Previous calculations also demonstrated that the stability of fully hydroxylated surface is 
comparable to the Al-terminated one and much more stable than O-terminated surface 28|, 



29J]. Similar to Al-terminated surface, graphene layer adsorbed on hydroxylated surface also 



belongs to vdW interaction. The most stable configuration is Bo (one surface O directly 
below the center of C-C bridge site), as shown in Fig. lc. The interlayer distance between 
graphene layer and hydroxylated surface ranges from 2.31 A to 3.18 A due to the rippled 
structure at the interface and the binding energy is -49.2 meV/C-atom, slightly smaller than 
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FIG. 2: HSE functional calculated electronic band structures of (a) single-layer graphene and (b) 
bilayer graphene on Al-terminated (0001) AI2O3 surface. The Fermi level is set to zero. The partial 
charge densities of VBM and CBM for the two structures are plotted as bottom and upper insets 
in (a) and (b), respectively. 



the binding strength of interlayers in graphite. The binding energy of Bo case is slightly 
lower than that of To and Ho within 1 meV/C-atom, indicating that these configurations 
are almost degenerate in energy. 

After knowing the stable interface structures, we turn to understand the effects of AI2O3 
substrate on the electronic property of graphene. The HSE band structure of single-layer 
graphene adsorbed on Al-terminated surface (Tai configuration) is shown in Fig. 2a. No- 
tably, the degeneracy of Dirac cone of graphene is lifted and a band gap of 182 meV (69 
meV in DFT+D/PBE calculations) appears at the Fermi level, generating an effective mass 
of the Dirac fermions ~ 8 x 1CT 3 m e . This is because the inversion symmetry in the 
graphene plane is broken on AI2O3 surface so the A, B sublattices, which is equivalent in 
free-standing graphene, is no longer equivalent. We find that the electrostatic potential 
difference of graphene C atoms at the A and B site differ by ~ 165 meV, close to the value 
of the band gap. The charge density distributions of conduction band minimum (CBM) 
and valence band maximum (VBM) (inset of Fig. 2a) are found separately located in the 
different A and B sublattices of graphene layer, which is consistent with our explanation. 

The sizes of band gaps of graphene in other high symmetrical configurations are close to 
that of Tai, indicating that a similar ON/OFF ratio could be obtained even if graphene slides 
away from its ground state Tai configuration. Remarkably, the band gap of graphene on Al- 
terminated AI2O3 is comparable to 10-nm-width GNRs in experiments 3 ij]. Meanwhile, the 
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effective mass of carriers in graphene/ A1 2 C>3 is an order lower than that of GNRs (~ 0.07 
m e ) from theoretical calculation^), strongly indicating the much higher carrier mobility 
of graphene/Al 2 03 in practice. The band alignment between the graphene and substrate 
surface is found to be type-I and there are no interface states around Fermi level for the Al- 
terminated interface (Fig. 2a). By distinguishing the respective graphene and AI2O3 states, 
we can calculate the band offset for this hybrid system. The valence band offset between 
graphene and AI2O3 is about 2.35 eV from the hybrid functional calculation (DFT+D/PBE 
level calculation seriously underestimates the valence band offset by 1 eV), comparable to 
that of silicon and AkO^flQl]. which is high enough for injection barrier. 

The electronic properties of bilayer graphene on Al-terminated surface, as shown in Fig. 
2b. The equilibrium graphene interlayer distance with AB stacking is 3.18 A. Similar to 
single-layer graphene, a band gap of 138 meV is induced in bilayer graphene due to the 
inversion symmetry breaking. The VBM and CBM are contributed by Ai and B 2 sublattices 
(Ai and Bi belong to the bottom layer while A 2 and B 2 belong to the top layer) in different 
graphene layers, respectively, as shown in Fig. 2b. The valence band offset between bilayer 
graphene and A1 2 C>3 is 2.26 eV, slightly smaller than that of single-layer graphene. It is 
important to notice that the band gaps of single-layer or bilayer graphene on Al-terminated 
surface are large enough to create sufficient ON/OFF ratio in FET applications {2 -5]. As 
we expected, the band gap of graphene on A1 2 C>3 is significantly larger than that on other 
popular substrates such as Si0 2 , Hf0 2 , and BN, indicating that a-Al 2 03 is a better choice 
for gate material 15 1 . 

The electronic properties of graphene on O-terminated surface are completely different 
from that of Al-terminated one. Due to the strong interaction between graphene and surface 
O atoms, the linear band characteristic of graphene disappears, as shown in Fig. 3a. The gap 
states in Fig. 3a originate from the hybridization between graphene C orbitals and surface O 
orbitals, shown as the inset in Fig. 3a. When a second graphene layer is placed on the surface 
(the graphene interlayer distance varies from 3.09 A to 3.75 A due to the rippled structure 
of the bottom layer), these gap states couple with the second layer graphene states and 
perturb its linear band distribution, as shown in Fig. 3b. A band gap of 140 meV appears 
for the second graphene layer, but the interface states exist inside the gap. Because these 
localized interface states can largely suppress the high mobility of graphene, O-terminated 
surface must be avoided to contact with graphene when gating graphene in experiments. 
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FIG. 3: HSE functional calculated electronic band structures of (a) single-layer graphene and (b) 
bilayer graphene on O-terminated (0001) AI2O3 surface; the partial charge densities of the partially 
occupied band at the Fermi level in (a) are plotted as inset. HSE electronic band structures of (c) 
single-layer graphene and (d) bilayer graphene on fully hydroxylated (0001) AI2O3 surface are also 
plotted. The Fermi level is set to zero. 



The electronic properties of graphene adsorbed on hydroxylated surface show similar 
behaviors to the Al-terminated surface, as shown in Figs. 3c and 3d for Bo configuration. 
No interface states around Fermi level are found and a band gap of 84 meV is induced 
in graphene layer with the effective mass of ~ 4 x 10~ 3 m e for the Dirac fermion. The 
valence band offset between graphene and hydroxylated surface is 3.39 eV, larger than that 
of Al-terminated one. Moreover, the electronic properties of To and Ho configurations are 
quite similar to that of Bo- Comparing to the single layer case, the bilayer graphene on 
the hydroxylated surface has a larger band gap of 126 meV and smaller band offset of 3.24 
eV, as shown in Fig. 3d. Clearly, O-terminated surface should be hydroxylated in order to 
make clean interface between graphene and AI2O3. Although the amorphous structure of 
Al 2 03 has been selected as a top gate in some experiments [8], it is difficult for us to directly 



7 



0.30 



^0. 



1.24 



& 0.18 



a 

cs 

CJ 

1 

3 

*5 0.12 



0.06 



a 




— ■- band gap 
— •— valence banc 


offset 











































-0.6 -0.4 41.2 0.0 0.2 0.4 
Electric field (V/A) 



2.6 



2.4 P 
0J 



2.2 



O 

T3 
C 



0.6 



1.8 



0.45 



O 0.36 



&0.27 



2.0 E "° 0.18 



0.09 



u 

D 






-•- band gap 
-•- valence band 


offset 
































■ 







-0.6 -0.4 -0.2 0.0 0.2 0.4 
Electric field (V/A) 



2.6 



2.4 £? 



2.2 SO 
O 



2.0 



73 
S 

1 



0.6 



1.8 



FIG. 4: HSE functional calculated band gap and valence band offset of (a) single-layer graphene 
and (b) bilayer graphene on Al-terminated AI2O3 surface as a function of external electric field 
(E ex t). The principal scheme of the computational model is shown as inset of (a). An E ex t is 
oriented normally to a surface and is assumed to be positive when it is directed upward. 

compare our results with these experiments due to the different phases. 

Finally, we considered the electronic properties of graphene/Al 2 03 hybrid system under 
an external electric field (E ext ) to simulate the gating effect in experiments 24| . Taking 
graphene on Al-terminated surface as examples, Fig. 4 shows the band gap and valence 
band offset of graphene/A^Os system as a function of E ext . For single-layer graphene on 
A1 2 3 , a negative E ext increases the band gap while the trend is opposite for a positive 
Eext, as shown in Fig. 4a. This is because a positive (negative) E ext increases (decreases) 
the distance between graphene and surface, e.g., the distance increases (decreases) from 
2.74 A to 2.86 (2.69) A when E ext = 0.6 (-0.6) V/A. The increase of the distance between 
graphene and substrate reduces the potential difference Aab, thus changing the band gap. 
The variation of potential difference also influences the valence band offset, as shown in Fig. 
4. The valence band offset increases (decreases) with a negative (positive) E ext and varies 
from 2.49 eV to 2.05 eV under -0.6 eV/A< E ext < 0.6 eV/A. 

The situation is quite different for bilayer graphene/Al 2 03 (Fig. 4b) under an E ext . Due 
to the symmetry, the band gap of bilayer graphene is determined mainly by the absolute 
difference of IAajbJ. 32| This value reaches minimum at E ext ~ 0.1 V/A when the system 
has a minimum anticrossing band gap around 112 meV. As E ext increases or decreases from 
this value, jA^Bal increases, which makes the band gap increases as well, as shown in 
Fig. 4b. For example, the band gap of bilayer graphene increases to 361 meV at E ext = 



8 



-0.4 V/A. The valence band offset changes between 2.37 eV and 1.93 eV for -0.6 eV/A< 
E ext < 0.6 eV/A. Thus, even under high E ext , the valence band offset (injection barrier) of 
graphene/Al 2 03 is still large enough. 

with physical insights and by first-principles calculations we demonstrate that a-A^O^ 
is an ideal high- A; gate material for graphene FETs. There are no interface states between 
graphene and Al-terminated or hydroxylated A1 2 3 surface and the valence band offset for 
this hybrid system is larger than 2 eV. Remarkably a band gap of ~ 180 meV appears 
in graphene layer deposited on Al-terminated surface, which is significantly larger than on 
other popular used substrates. Furthermore, the band gap of graphene/Al 2 0,3 system could 
be tuned by an external electric field for practice applications. 

The authors thank the helpful discussion with Prof. Jaejun Yu (SNU). The work at 
NREL was supported by the U.S. Department of Energy under Contract No. DE-AC36- 
08GO28308. 
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